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A theory is developed for interband tunneling in semiconducting carbon nanotube and graphene 
nanoribbon p-n junction diodes. Characteristic length and energy scales that dictate the tunneling 
probabilities and currents are evaluated. By comparing the Zener tunneling processes in these struc- 
tures to traditional group IV and III-V semiconductors, it is proved that for identical bandgaps, 
carbon based ID structures have higher tunneling probabilities. The high tunneling current magni- 
tudes for ID carbon structures suggest the distinct feasibility of high-performance tunneling-based 
field-effect transistors. 

PACS numbers: 73.40.-c, 73.40.Kp, 73.40.Lq 

Keywords: interband, Zener, tunneling, graphene, nanoribbon, carbon, nanotube, diode, transistor 



Carbon-based ID materials such as nanotubes (CNTs) 
and graphene nanoribbons (GNRs) are currently under 
extensive investigation for the novel fundamental physics 
they exhibit, as well as possible applications they might 
have in the future^. A large class of traditional semi- 
conductor devices rely on the quantum mechanical tun- 
neling of carries through classically forbidden barriers. 
Among these, the Esaki diode, resonant tunneling diode, 
and backward diode are the prime examples^. In ad- 
dition, the high-field electrical breakdown in a number 
of semiconductors occurs by interband Zener tunneling. 
The phenomena of tunneling has been studied extensively 
for traditional parabolic-bandgap semiconductors in 3D 
bulk, as well as quasi-2D and quasi- ID heterostructures. 
Semiconducting CNTs and GNRs do not have parabolic 
bandstructures, and the carrier transport in them ap- 
proaches the ideal ID case. In that light, it is timely to 
examine the phenomena of tunneling in these materials. 

Tunneling currents in semiconducting CNT p-n 
jun ction s have been measured and analyzed recently 
(se d 4 * 5 * 6 j). For tunneling probabilities, an energy- 
dependent carrier effective mass has been used earlier^, to 
take advantage of previously existing results of parabolic 
bandstructure semiconductors. In this work, we evalu- 
ate the Zener tunneling probabilities of CNT and GNR 
based p-n diodes starting from their intrinsic bandstruc- 
tures, which removes the need to define an effective mass. 
In addition to interband tunneling probabilities, a num- 
ber of fundamental associated parameters characterizing 
the tunneling process are found. 

The bandstructure of the nth subband of a semicon- 
ducting CNT or GNR is given bj0 

£ = shvF\fk1 + k%, (1) 

where 2irh is the Planck's constant, and Vp ~ 10 8 cm/s 
is the Fermi velocity characterizing the bandstructure of 
graphene. s = +1 denotes the conduction band, and s = 
— 1 denotes the valence band. The electron momentum 
along the CNT or GNR axis is hk x . 

For GNRs, the transverse momentum is quantized 



by the ribbon widtfP; k n = nir/3W where n = 
±1, ±2, ±4, ±5, ±7, ±8... for a GNR of dimensions 
(x,y) = (L, W) where W << L. The corresponding 
bandgap is £ g — 2Uvpk\ — 2nhvp/3W ~ 1.3 j"W eV 
(where W is in nm). For comparison, a semiconducting 
CNT of diameter D has the same bandstructure with 
fci = 2/3D, and a bandgap of £ g = 4hv F /3D ~ 0.8/ D 
eV, where D is in nm. If W = ttD/2, the properties 
(bandgap, bandstructure) of semiconducting CNTs and 
GNRs are similar. The results derived below are applica- 
ble to GNRs and CNTs on equal footing. It is assumed 
that the length of the GNR (CNT) is much larger than 
the width (diameter) such that the longitudinal momen- 
tum of carriers in the ribbon are quasi-continuous. 

We now evaluate the interband tunneling probability 
inan+- p+ GNR or CNT diode of bandgap £ g . We 
consider that the doping of the n and p-sides are such 
that the equilibrium Fermi level is at the conduction band 
edge (£ c ) in the n-side, and at the valence band edge (£ v ) 
on the p-side. Such doping could be either chemical, or 
electrostatic^. Under this situation, a forward bias would 
not lead to current flow (this is similar to the 'backward 
diode 251 ). When a reverse bias eV is applied, the band 
diagram looks as shown in Figure [T] Let the electric 
field in the junction region be T. We assume that the 
depletion region thickness, and the net electric field does 
not change appreciably from the equilibrium values (true 
under small bias voltages). Then the potential energy 
barrier seen by electrons in the valence band of the p + 
side is 

V (x) = eTx, (2) 

in the range < x < d, such that d = £ g /eT . This 
is indicated in Figure [T] £ a is the Dirac point of the 
underlying graphene bandstructure, for both GNRs and 
CNTs, and serves as a convenient reference of energy in 
the problem. 

Since the energy of carriers near the band edge (k x 0) 
is conserved during the tunneling of electrons from the 
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FIG. 1: Interband tunneling in a reverse-biased GNR/CNT 
p-n junction, and the potential barrier seen by tunneling elec- 
trons. 



valence to the conduction band, the condition 
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holds for the wavevector k x at all x. Within the tunneling 
barrier, the wavevector is imaginary. Denoting this by 
k x = in x where n x is real, we obtain 



k x (x) = k Q 



(4) 



where fco = £ g /2hvF is a characteristic wavevector for 
tunneling. Since J* y/1 - (1 - |) 2 = 7rd/4, the WKB 
band-to-band tunneling (BTBT) probability for the p + — 



junction given by T ~ exp [— 2| J Q +d n x (x)dx\] leads to 
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which can be expressed as IV^^ ~ exp(— J-'q/J 7 ), 
where J-b = ff£g I '^ehvp denotes the characteristic elec- 
tric field at the junction for the onset of strong tunneling. 
The corresponding characteristic barrier thickness for the 
onset of strong tunneling is do ~ £ g jeT§ — 4hvp /Tr£ g . 

Using the value of the Fermi velocity, the characteristic 
field evaluates to J-q ~ 12.6 x (£ g ) 2 MV/cm, and the char- 
acteristic tunneling distance evaluates to e?o ~ 0.8/ £ g nm, 
where £ g is the bandgap of the GNR or CNT expressed 
in eV in both these expressions. As an example, for a 
GNR with W = 5 nm, £ g ~ 0.275 eV, the characteristic 
tunneling field is ~ 0.9 MV/cm, and the characteristic 
tunneling barrier thickness is ~ 11 nm. 

The tunneling probabilities evaluated above are for the 
first subband (n = 1). For the nth conduction and va- 
lence subbands of CNTs and GNRs, the effective subband 



gap scales as ~ n£ g . The tunneling probabilities of car- 
riers near the respective subband-edges is then given by 
Twkb,u ~ exp(— irhvpkn/ J 7 ), which decay as exp(— n 2 ), 
indicating a rather strong damping the tunneling proba- 
bilities of higher subbands. This result turns out to be 
identical to one for the nth transverse mode of a zero-gap 
2D graphene p-n junction, as first evaluated by Chcianov 
and Fal'ko (see^). Except for the narrowest bandgap 
CNTs and GNRs, the tunneling is primarily from the 
1st valence subband in the p-side to the 1st conduction 
subband on the n-side. 

Note that the above derivation uses a triangular bar- 
rier approximation. It has been shown by Kane^ 2 that a 
parabolic barrier more accurately represents the physics 
of the tunneling process, and leads to an exponential fac- 
tor with a different coefficient than from the triangular 
barrier approximation. The difference is small, for the 
rest of this work, we will use the result derived above. 

The tunneling current may now be written in the form 
I T = e(g s g v /L)J2 k Jfv-fc\TwKB'xv g (k x ), where g s = 2 
is the spin degeneracy, g v is the valley degeneracy (= 2 
for CNTs and = 1 for GNRsP), L is the length of the 
CNT or GNR, f v ,f c are the occupation functions of the 
valence band states in the p-type side and the conduction 
band in the n-side respectively, and v g (k x ) is the group 
velocity. This sum may be converted into the equivalent 
integral 
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where f v (£) = 1/(1 + exp [(£ - eV)/k B T}) and f c (£) = 
l/(l + exp [S/ksT]). Here ks is the Boltzmann constant. 
For the band diagram shown in Figure [l] the net tunnel- 
ing current evaluates to 

It =^-T WKB xV T ln[hl + cosh ^)], (7) 

where Vt — ksT/e. This expression captures the tem- 
perature and bias voltage dependence of the tunnel- 
ing current. If the applied bias is much greater than 
the thermal energy (V » Vt), the current reduces 
to a form similar to the Landauer expression It ~ 
(2g v e 2 /h)T WKB (V - Vyln4). Thus, the tunneling cur- 
rent has a negative temperature coefficient at high bias 
conditions. The dependence of tunneling currents in 
GNRs on voltage and temperature are illustrated in Fig- 
ure |2| (a, b). 

The tunneling current per unit width of GNRs is max- 
imized for Wq = y /, 2n 3 hvp/9eJ-. For example, for T ~ 1 
MV/cm, the GNR width for maximum current drive is 
Wq ~ 6.5 nm, and the current density at that width 
is ~ 450 /iA//im, comparable to traditional FETs. For 
thinner ribbons at higher fields, the current densities can 
approach 1000 /iA//im. The dependence of the current 
densities on temperature, electric field, and GNR widths 
are shown in Figure 2](c, d). 

The BTBT probability in traditional parabolic- 
bandstructure semiconductors in the triangular barrier 
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FIG. 2: Tunneling currents in GNR p-n junctions: (a) W =5 
nm device, voltage dependence at various temperatures, (b) 
temperature dependence at various voltages. Tunneling cur- 
rents per unit width for different GNR widths at (c) various 
temperatures at T = 1 MV/cm, and (d) for various T at 300 
K. To maximize the tunneling current density, an optimum 
(J- ', W) combination exists, as illustrated by the dashed lines. 



approximation depends on the bandgap as^ T V araboiic ~ 

exp (-AV2m*£g /2 /SetLF), where m* is a reduced carrier 
effective mass. The BTBT probability for GNRs and 
CNTs retains the same dependence on electric field, but 
due to the difference in bandstructure there is a stronger 
dependence on the bandgap. The effective mass does ap- 
pear in the tunneling probability of the GNR or CNT 
diodes since their bandstructure is not parabolic at any 
energy. If one compares the Zener tunneling probabili- 
ties in diodes made of CNTs or GNRs with other direct- 
bandgap semiconductors of the same bandgap, then the 
ratio 

Tcarbon r / 7r \/^9 4i/2m* . -. , . 

exp [- -r-, -(— — ) (8) 

J- parabolic enJ- Av F 3 

indicates that the GNR or CNT p-n diode will have 



a higher interband tunneling probability if the relation 
£ g < (16t>i?/37r) 2 x2m* is satisfied. From the k ■ p theory 
for traditional direct-bandgap semiconductors, the elec- 
tron effective mass (in the conduction band) is related to 
the bandgap by the approximation^ to* w (£ g /20)mo, 
where the bandgap is in eV and too is the free electron 
mass. This leads to the requirement 

m v 2 F >A 2 xlOeV, (9) 
16 

which is always satisfied since the LHS is n%oVp ~ 5.7 eV 
and the RHS is 3.5 eV. Thus, the CNT or GNR p-n diode 
will always have a higher reverse-bias Zener tunneling 
probability than a traditional semiconductor of the same 
bandgap. 

Two more facts tilt the tunneling probability decisively 
in favor of CNTs and GNRs. First, for bulk 3D p-n 
junctions, the transverse kinetic energy of carriers can 
be large, and leads to a further exponential decrease of 
carrier tunneling probability 3 , which is avoided in ID 
structures. Second, if normal parabolic bandgap semi- 
conductors are shrunk to length scales comparable to 
that of CNTs and GNRs, their bandgap and the effective 
masses increase further due to quantum confinement. 

Though Zener tunneling currents are detrimental in 
traditional devices such as rectifiers, field-effect and bipo- 
lar transistors, it is important to note that the funda- 
mental switching action in these devices is controlled by 
thermionic emission over barriers, which require a min- 
imum of (fc^T/e) ln(10) ~ 60 mV per decade change of 
current at 300 K (the 'classical' limit). However, a new 
crop of tunncling-FETs have been recently proposed and 
demonstrated 4 ^ 14 ! 15 '), which rely on the very mechanism 
studied in this work. These devices are capable of reach- 
ing far below the classical limit for switching by exploit- 
ing quantum mechanical tunneling, ft is for such devices 
that high interband tunneling current drives in carbon- 
based ID nanostructures hold a distinct advantage, and 
much promise in the future. 
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